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ABSTRACT 

The general properties of the non-thermal non-pulsed X-ray emission of 
rotation powered pulsars are investigated in the context of a pulsar wind nebula 
model. An examination of the observed X-ray emission from a sample of 23 
pulsars in the energy range between 2-10 keV reveals that the relation of X-ray 
luminosity, L^., to the pulsar spin down power, E, is steeper for the non-pulsed 
component than for the pulsed component. Specifically, L"^"' oc £'i-4±o.i 
the non-pulsed component, whereas L^"' oc iji-2±008 f^^. ^^le pulsed component. 
The former relation is consistent with emission from a pulsar wind nebula 
model in which L"^"' cx E^^'^ where p is the power law index of the electron 
energy distribution. The relation for the pulsed component, on the other hand, 
is consistent with a magnetospheric emission model. In addition, the photon 
spectral index, F, was found to be correlated to the conversion efficiency of 
spin down power to non-pulsed X-ray emission with greater efficiencies for 
r ~ 2 — 2.5 than for F ~ 1.5 — 2. Such a relation is naturally understood 
within the framework of a pulsar wind nebula model with the former relation 
corresponding to the emission of X-rays in the fast cooling regime and the latter 
relation corresponding to emission in the slow cooling regime. 

The X-ray properties of pulsar wind nebulae are sensitive to the physical 
conditions (e.g., the density and magnetic field) of the interstellar medium and 
can lead to important differences between the X-ray emission characteristics 
(luminosity, photon spectral index and emission morphology) of pulsars in 
various environments. Such wind nebulae can contribute to the non-thermal 
symmetric emission morphology (point-like) and elongated emission morphology 
(tail-like) from sources similar to Geminga and PSR B1757-24. 

Subject headings: stars: neutron - X-rays: stars - pulsars: general - radiation 
mechanisms: non-thermal 
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1. INTRODUCTION 

The study of the emission characteristics from rotation powered pulsars has been 
a subject of long standing interest. Our knowledge of the fundamental properties (e.g., 
mass, spin, and magnetic field) of the underlying neutron star stems from detailed spectral 
and timing investigations. To facilitate an understanding of the mechanism by which the 
loss of rotational energy is converted into high energy radiation, many observational and 
theoretical studies have sought to determine the relationship between the X-ray luminosity, 
Lx, and the rate of rotational energy loss or spin down power, E. Indeed, a correlation of 
the form Ly; oc E^-^^ was found in Einstein data by Seward & Wang (1988). Subsequent 
studies using a larger sample of pulsars led to a relation of the form, (x E based on 
ROSAT data (see Becker & Triimper 1997) and oc E^'^ based on ASCA data (see Saito 
1998). Recently, a reanalysis of 39 pulsars based on data obtained from several X-ray 
satellites by Possenti et al. (2002) led to an intermediate relation oc similar to 

that of Seward & Wang (1988). 

While the deduced existence of a correlation between and E suggests that the 
observed X-rays are produced by a process which taps the rotational energy of the neutron 
star, a detailed description of the mechanism remains elusive. This is, in part, a result 
of the fact that the data from different sateUites are obtained in different energy ranges. 
The results can be affected, for example, by interstellar absorption especially for those 
pulsars studied in the soft X-ray regime of ROSAT (0.1- 2.4 keV). In addition, the total 
X-ray luminosity is composed of contributions from both the pulsed and non-pulsed 
components, and these components are likely to reflect physical conditions in diverse spatial 
environments. 

The X-rays radiated by rotation powered pulsars could include contributions from 5 
distinct components with different origins: 

1. Non-thermal non-pulsed diffuse radiation can be emitted from a pulsar wind nebula. 
The radiation will be discussed in §2, and we suggest that this emission component is the 
major contribution to the non-pulsed X-ray radiation in the ASCA observations of pulsars. 

2. Non-thermal non-pulsed radiation from the pulsar magnetosphere could be important as 
suggested by Becker et al. (2004) if the angular resolution of X-ray detector is sufficiently 
high to separate this contribution from that of the nebula. Although such a component may 
be present in the ASCA data that we use, due to its low angular resolution, it is not hkely 
to dominate the contribution from the pulsar wind nebula. For example, Tennant et al. 
(2001) have found a non-thermal non-pulsed component from the Crab pulsar by using the 
High-Resolution Camera of the Chandra X-ray Observatory. They estimate that the X-ray 
luminosity of this component is ~ 10^"^ergs~^, whereas, the non-thermal non-pulsed X-ray 
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luminosity from the Crab nebula detected by ASCA (Saito et al. 1997a) is ~ lO^^ergs"^. 
The specific origin of the non-thermal non-pulsed component, not associated with the 
nebula, is not clear. Although it contaminates the emission of the nebula, it is small. 
Provided that this contribution is small for all the pulsars we study in this investigation, 
this component should not affect our analysis. 

3. The non-thermal pulsed radiation component is generally believed to be produced from 
the pulsar magnetosphere. This follows from the fact that the motion of charged particles 
within the light cylinder is strongly affected by the magnetic field, and the radiation 

of charged particles in the open field hues should be pulsed. This component could be 
produced in the vicinity of the polar cap as a result of inverse Compton scattering of higher 
order generation pairs of particles on soft photons emitted by the neutron star (Zhang 
& Harding 2000) or in the outer magnetosphere as a result of synchrotron radiation of 
downward cascades from the outer gap electron/positron particles (Cheng & Zhang 1999). 

4. The thermal non-pulscd radiation from the pulsar surface can contribute to the soft 
X-ray bands, typically, characterized by kT ~ 0.1 keV (Cheng & Zhang 1999). In this 
paper, we have used the ASCA data on pulsars so this component will not contribute to the 
2-10 keV ASCA energy band. 

5. The thermal pulsed radiation from the pulsar surface likely originates in a polar cap 
on the neutron star surface and can contribute to the hard X-ray band {kT ~ 1 keV). In 
analyzing the pulsed emission from pulsars, we follow the model of Cheng & Zhang (1999) 
where the thermal pulsed X-ray emission is included in their model prediction. 

Thus, the X-ray emission observed from pulsars in the ASCA sample primarily consists 
of contributions from components 1, 3, and 5. We have assumed that component 2 is 
not important and component 4 docs not significantly contribute to the 2-10 keV energy 
band. On the other hand, Becker & Triimper(1997) have used ROSAT data and obtained 
a linear relation between the X-ray luminosity consisting of components 2 and 3 and 
spin-down power. In the present paper, we concentrate on the X-ray properties in the 
energy band of ASCA (2-10 keV). The X-rays, then, consists of a thermal component from 
the polar cap, a non-thermal pulsed component from the magnetosphere and a non-thermal, 
non-pulsed emission from the nebula. That is, we have explicitly assumed that the pulsar 
magnetosphere and the pulsar wind nebula are primarily responsible for the pulsed and 
non-thermal non-pulsed component respectively. 

Since the pulsed and non-pulsed emission have different origins, the relationship 
between the total X-ray luminosity and spin down power is not expected, in general, to be 
represented by a single power law. In fact, the non-thermal pulsed and non-pulscd emission 
components are distinguished by different spectral signatures and conversion efficiencies. To 
determine their relative importance and their relationship to the spin down power, we have 
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reexamined the pulsar data obtained in the 2-10 keV energy band of the ASCA sateUite, 
focusing on separating the non-pulsed and pulsed emission data. In the next section, we 
provide a simplified description of the emission from a pulsar wind nebula, based upon 
earlier work by Chevalier (2000), where the relativistic pulsar wind interacts with the 
interstellar medium. In §3, the X-ray properties of the pulsed and non-pulsed components 
are collected and the relation between their luminosities with spin down power and photon 
spectral index are determined. With these results in hand, the pulsar wind nebula model is 
used to explain the occurrence of the tail-like emission morphology in Geminga and PSR B 
1757-24 in §4. Finally, we summarize and conclude in the last section. 



2. THEORETICAL FRAMEWORK 

A theoretical description of the interaction between a pulsar and its nebula, as applied 
to the Crab nebula, was first outlined in a seminal paper by Rees & Gunn (1974). In 
this model, the central pulsar generates a highly relativistic particle dominated wind 
which passes through the medium in the supernova remnant, forming a shock front. The 
electrons and positrons in the shock are envisioned to be accelerated to a power law energy 
distribution and to radiate synchrotron photons in the downstream region. 

Since the magnetic field in the nebula derives from the pulsar, a magnetohydrodynamical 
(MHD) model, rather than a hydrodynamical model, of the pulsar wind nebula was required. 

Siich a model, in the steady state and spherically symmetric approximation, was developed 
by Kennel & Coroniti (1984) for the Crab nebula. A hot relativistic positronic plasma flow 
was assumed to be terminated by a strong MHD shock, decelerating the flow and producing 
a non-thermal distribution of electrons. The magnetization parameter, cr, introduced to 
describe the efficiency for conversion of energy contained in the pulsar wind into synchrotron 
radiation, was represented by the ratio of the electromagnetic energy flux to the particle 
kinetic energy flux deflned as: 

Here n is the particle number density, B is the magnetic held, 7^ is the Lorentz factor of 
relativistic particles in the wind, m is the particle mass, and c is the speed of light. As 
pointed out by Kennel & Coroniti (1984), the magnetization parameter must be small (for 
the Crab nebula, a ~ 0.003) in order that sufficient compression occurs in the shock for 
transformation of the bulk fiow energy to random particle motion and to the subsequent 
production of the observed synchrotron radiation. Since the magnetic energy density must 
dominate the particle energy density just inside the light cylinder, the wind nebula must 
have evolved from a high cr to a low a state to produce the observed properties. The 



- 5 - 



manner in which this transformation takes place is central for detailed pulsar wind models 
and is not understood in detail. Nevertheless, the low a model can account for the basic 
high energy properties of the Crab, including the position of the optical wisps, the spectrum 
of the nebula, and the size of the nebula as a function of wavelength (van den Bergh & 
Pritchet 1989). 

In the present paper, we concentrate on the X-ray properties of the pulsar wind 
nebula, using a simple one zone model similar to that developed by Chevalier (2000) 
for a description of the wind interaction with the surrounding medium in a variety of 
environments. Although the model does not provide a description of the spatial details 
associated with equatorial tori and polar jets observed in Chandra studies of young 
pulsar wind nebulae, it can explain the global X-ray properties of the nebula to a first 
approximation. 



2.1. One-zone model of X-rays from wind nebulae 

In general, the energy in the shock waves is stored in the magnetic field as well as 
in the proton (ion) and electron particle components. We assume that the fractional 
energy density of the magnetic field, e^, is ~ 10~^ — 10^^ (Kennel & Coroniti 1984). 
Provided that the spin down power is eventually carried away by pulsar wind particles and 
assuming equipartition of energy between the electrons and protons, the fractional energy 
densities in the proton and electron components should satisfy ~ eg 0.5. For a given 

in the shock, the magnetic field at the shock termination radius, Rg is estimated as 
B = {QesE / Rlcf'^ , where E = AnlP/P^ is the spin down power of the pulsar. Here, P 
and P are the spin period and its derivative, and / is the moment of inertia of the neutron 
star (~ lO'^^g cm^). 

At the shock front, the electrons attain a power law distribution corresponding to 
A^(7) oc for 7m < 7 < 7maxj where 7 is the Lorentz factor and 7^ = '^^ec^/u;- An 
estimate for 7max can be obtained by equating the synchrotron cooling timescale to the 
electron acceleration timescale. The former timescale is given by tgyn = Qirmec/aTjB^, and 
the latter timescale is given by tacc = 'jniec/eB, leading to 7max = (Ovre/crri?)^/^, where nie 
is the electron mass and ar is the Thompson cross section. The synchrotron power of an 
electron with Lorentz factor 7 is (Blumenthal & Gould 1970) 

P(7) = larc^'^. (2) 
Following Chevalier (2000), the total rate of particles ejected by pulsars at the Lorentz 
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factor 7 is (Chevalier 2000) 

iV(7) = {p- l)Y^-\^^mcY'Er''. (3) 

The resulting total number of radiating particles becomes A'"(7) ~ N{'y)t, where i is a 
characteristic timescale. For a young pulsar with a supernova remnant similar to the Crab, 
t is identified with the age of the nebula while for old pulsars or millisecond pulsars, t is 
comparable to the flow timescale within the nebula. Thus, at the emitting frequency of 
an electron 2/(7) = '-j'^eB /2Txm^c, the luminosity of the particle is about P{'~f)N{'~f) which 
should be less than E. The detailed properties of the X-ray nebula are a function of Cg, e^, 
Rs, E, and the two critical frequencies z/^ and Uc where is the frequency radiated by the 
electrons with the Lorentz factor of 7^, and 

i/c = ISneniec/aTt'^B^ (4) 

is the electron synchrotron cooling frequency (see Chevalier 2000 and references therein). 

The X-ray luminosity and spectral index depend on the inequality between Ux and 
Uc- For typical values of ~ 0.5, 7^ ~ 10^ and B ~ 10^'^ G, z/^ 10^^ Hz is always less 
than i>x- In this case, the X-ray luminosity from the total emission of particles per unit 
frequency can be calculated according to the synchrotron spectral profiles which have been 
presented in Figure 1 of the paper by Sari, Piran & Narayan (1998). For ux > i^c, the 
luminosity per unit frequency can be obtained following equations (7) and (8) of Sari, Piran 
& Narayan (1998) 

~ {uJu,)-'/\u/uJ-P/^L,,^^, (5) 

where -Li/,max = -^(7m)-P(7m), and u ~ 10^^ Hz. This expression corresponds to the fast 
coohng regime, and is equivalent to equation (9) in Chevalier (2000) reflecting our use of 
the same form of particle number injection rate by a pulsar (see equation 3 in this section). 
On the other hand, if vx < (denoted as slow cooling), the synchrotron X-ray luminosity 
is estimated as (Sari, Piran & Narayan 1998) 

~ {iy/iym)-^-^^^^L„. (6) 

It can be seen that the X-ray luminosity in the fast cooling regime is larger than that in 
the slow cooling regime by a factor (I'c/i^xY^'^ if i^x is in the slow cooling regime. Since 
the range of p due to shock acceleration lies between 2 and 3 (sec Achtcrbcrg ct al. 2001; 
Lemoine & Pelletier 2003 and references therein), the photon index in the slow cooling 
regime is expected to lie in the range from 1.5 to 2 whereas the photon index in the fast 
coohng regime is expected to be in the range of 2 to 2.5. 
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2.2. Interactions between the pulsar wind and surrounding medium 

Within the theoretical framework, the energy of the electrons is converted to X-ray 
radiation solely at the shock termination radius. For pulsar motion which is subsonic, the 
determination of this radius is obtained via the balance between the wind ram pressure and 
total magnetic and particle pressure within the nebula (Rees & Gunn 1974). In this case, 
the wind bubble will be nearly centered about the position of the pulsar. Rees & Gunn 
(1974) estimated i?^ ~ 3 x 10^^ cm, which is consistent with the size of the inner X-ray ring 
of the Crab nebula (~ 0.1 pc, see Kennel & Coroniti 1984; Weisskopf et al. 2000). The 
termination shock picture is not necessarily restricted to slowly moving pulsars since even 
for pulsars moving at several hundred km s^^ the motion can be subsonic in, for example, 
regions where the gas has been heated to temperatures of ~ 10® K by a supernova shock. 
The radius of the termination shock can be estimated as: 

Rs^{^f'-^x^^''E'J,'B;;^a^u.^ (7) 

where E:^,.l is the pulsar spin-down power in units of 10^^ ergs s~^ and -BmG is the magnetic 
field strength in the nebula in milligauss. With the observed values of the Crab pulsar and 
its nebula (£^34 = 5 x 10"^ and B^q — 0.5), consistency of the shock radius is easily achieved. 

For supersonic motion, the nebula will form a bow shock morphology. In this case, 
the termination shock radius is given by the balance of the ram pressure between the wind 
particles and the medium at the head of the bow shock: 

E 1 2 



2AsM^^p, (8) 



where pisM is the density of the interstellar medium and is the pulsar's proper motion 
velocity. The termination shock radius is given as: 

Rs ^ ( , ^ , ~ 3 X 10^«£;3f nr^/^^-JooCm, (9) 
STrpisMt'pC 

where ni, and t'p,ioo a-re the number density in the interstellar medium in units of 1 particle 
cm~^ and the pulsar space velocity in units of 100 km s~^. As an example, Caraveo et al. 
(2003) recently discovered a bow shock structure coincident with the Geminga pulsar based 
on observations obtained with XMM-Newton. An estimate of the termination radius for 
Geminga follows from its spin down power, £■ ~ 3.2 x lO^^erg s~^, distance of ~ 160 pc, 
and proper motion velocity of Vp ~ 120 km (Bignami & Caraveo 1993), leading to a 
radius of 4 x 10^^ cm which is consistent with the observational constraint on Geminga's 
compact X-ray nebula. 
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The general features of the emission region surrounding rotation powered pulsars in 
the Galaxy will be bracketed by these two wind nebula morphologies. Since the range in 
magnetic field strengths in the shock can vary from ~ 10~^ G to as high as 10~^ G in 
the Vela nebula and 5 x 10"^ G in the Crab nebula for the very young pulsars within 
supernova remnants, a range in the conversion efficiency of the spin down energy into 
nebula luminosity is expected. The relation between the nebula X-ray luminosity and the 
pulsar's spin down power will be deferred until §3. 

2.3. General properties of the wind nebulae in X-rays 

In the above we have discussed the pulsar wind nebula models in different contexts 
within a one- zone approximation. The luminosities and spectral properties of the nebulae 
can be determined by Eqs. (5) and (6) in the fast and slow cooling regime respectively. In 
the fast cooling case, the efficiency of converting spin down power to nebula luminosity is 
higher than that in the slow coohng case. Assuming that eg ~ 0.5, €b ~ 0.01, ^-^ ~ 10^ and 
p 2.2 (see Bednarz & Ostrowski 1998), the ratio of the X-ray luminosity, ~ lyL^,, to 
the spin down power at ~ 10^^ Hz can be estimated as Lj./E ~ 10~^ — 10~^ for the fast 
cooling case and L^/E ^ 10^^ — 10^"^ for the slow cooling case if t ~ 10^ s, and B ~ 10~^ 
G. The spectral properties in these two regimes differ since the photon index in X-rays is 
r = (p + 2)/2, and F = (p + l)/2 for fast and slow cooling respectively. Since p varies from 
2 to 3, r ~ 1.5 — 2.0 (slow cooling) and F ~ 2.0 — 2.5 (fast cooling). These are typical of 
the photon indices of observed X-ray nebulae (see Table 1 of Chevalier 2000). 

As described in §2.1, the fast and slow cooling regimes are distinguished by the X-ray 
and cooling frequencies. Since the cooling frequency is dependent on the magnetic field and 
the flow time scale of the nebula, it is possible that different regions of the nebula may 
correspond to different cooling regimes. For example, the inner region of the nebula may 
correspond to the slow cooling regime, whereas the fast cooling may be more appropriate 
in the outer region of the nebula. In this case, the power law index of the X-ray spectrum 
will steepen by ~ 0.5 between the inner and outer region. This may be consistent with the 
spectral results from the high spatial resolution observations of the nebulae surrounding 
PSR B1823-13 obtained with the XMM-Newton satellite (Gaensler et al. 2003a). 

Because the spatial resolution of current X-ray detectors is limited the emission from 
compact nebulae of size ~ Rg may be contaminated by the non-thermal X-ray emission 
from the pulsar magnetosphere. However, the pulsar magnetospheric radiation is assumed 
to be pulsed, whereas the nebula contribution is non-pulsed. Therefore, the separation of 
the non-thermal X-ray emission into its pulsed and non-pulsed components is essential for 
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deter mining their relative contributions. 

3. PULSED AND NON-PULSED EMISSION 

Since the pulsed and non-pulsed emission components originate from different regions, 
their X-ray luminosities are not expected to exhibit a similar relationship to the spin 
down power. For example, Wang & Zhao (2004) recently found that the correlation of the 
observed pulsed X-ray luminosity with spin down power is not well matched to the linear 
form found by Becker & Triimper (1997) for the total luminosity in the ROSAT energy 
band (0.1-2.4 keV, see Figure 8 of Wang & Zhao 2004). We note, however, that studies in 
this lower energy band are incomplete since pulsars can emit a significant fraction of their 
luminosity at higher energies. To examine the X-ray luminosity- spin down power relation 
at higher energies (2-10 keV), Kawai et al. (1998) selected a sample of bright nebulae and 
found a relation in the form Lx oc £'i-27±o.i7 gj^cg ^j^g individual components were not 
separately analyzed, the result was insufficient to establish the relation between the pulsed 
or non-pulsed emission with spin down power. 

3.1. Lx — E relations from ASCA data 

Over 50 rotation powered pulsars have been detected in the X-ray band (Becker & 
Aschenbach 2002), but only a fraction of them have been resolved with pulsed non-thermal 
components. Here, we have selected 23 X-ray pulsars with both pulsed and non-pulsed 
emission measurements obtained from the ASCA mission (hsted in Table 1). The pulsed 
and non-pulsed X-ray luminosities of these selected pulsars are directly taken from the cited 
references. This pulsar sample includes 19 normal pulsars and 4 millisecond pulsars whose 
observational data are taken from the recent references (see Table 1). Since the ASCA 
satellite does not have high spatial resolution, the X-ray luminosity of the pulsars in the 
ASCA field is composed of emission from the pulsar's magnetosphere and compact pulsar 
wind nebula. For reference, the total pulsed plus non-pulsed X-ray luminosity in the ASCA 
energy range (2-10 keV) is plotted versus spin down power in Figure 1. A correlation is 
found which is consistent with the form oc E^-^ found by Saito (1998), but the best fit 
form of this correlation is found to be oc £'i-35±o.2 fjgj-g^ ^j^g g^ror in the power law 
exponent represents ilcr corresponding to the scatter in the observed data points, which 
may reflect variations in ee, e^, 7w,, and uncertainties in distance from pulsar to pulsar. 
Our best fit power law relation is consistent with the conclusion of Possenti et al. (2002), 
who used a sample of 39 pulsars observed mainly by ROSAT and data from ASCA, RXTE, 
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BeppoSAX, Chandra and XMM-Newton. Although Possenti et al. (2002) considered their 
result (La; oc E^-^"^) statistically Tinacceptable, our result is not subject to the uncertainties 
associated with the normalization of different satellite data to obtain the X-ray luminosity 
between 2-fO keV (cf. Table f of Possenti et al. 2002), for which the extrapolation relied 
on the uncertain photon index in the ROSAT energy band. 

The X-ray luminosity associated with the pulsed emission component is illustrated 
versus spin down power in Figure 2. A correlation separate from the total luminosity is 
found which is inconsistent with either the form oc E or oc ^3/2. The best fitting 
function to the data is found to be Lx,pui — (1-0 ± 0.6) x iQ-^^E^-^^^^-O^^ which significantly 
deviates from the 3/2 power law relation proposed by Saito (1998). Such a relation is 
consistent with the relation oc E^-^^ derived from the theoretical X-ray magnetospheric 
emission model of Cheng & Zhang (1999), however, this latter result is not without 
uncertainties since the inclination angle of the magnetic field with respect to the rotation 
axis and the viewing angle are not well determined. The observed conversion efficiencies 
are found to range from ~ 10~^ — 9 x 10~^, which is not in conflict with model predictions 
(Cheng, Cn & Zhang 1998; Cheng & Zhang 1999). 

A correlation is also found to exist between the non-pulsed X-ray luminosity and 
the spin down power as shown in Figure 3. The data points are consistent with the 
previous ASCA relation of the type oc E^^"^, but the correlation is also consistent with 
-f'x.npui oc iJi-4±o.i^ Upon comparison to the results from Fig. 1, this power law relation is a 
consequence of the fact that the pulsar emission in the ASCA sample is dominated by the 
non-pulsed radiation component. The conversion efficiency for the non-pulsed component 
overlaps with that for the pulsed component, but extends to efficiencies as high as 0.1. 

Generally, the X-ray luminosity of pulsar wind nebulae (Chevalier 2000) is a nonlinear 
function of the spin down power. As can be seen from equation (9) of Chevalier (2000), the 
nebula's X-ray luminosity follows from 

Here, eg and €b are assumed to be constant, but Rg oc E^^'^ from Eqs. (7) and (9). Although 
the above equation corresponds to the fast cooling regime, the dependence of Lx on E for 
the slow cooling regime is unchanged because the ratio of the X-ray luminosities in these 
two regimes depends only on the cooling frequency z/g, which is independent of E. We note 
that 7^ also depends on the spin down power so the explicit dependence of on E remains 
to be theoretically determined. To estimate this dependence we make use of the results of 
Ruderman (1981) and Arons (1983), who argued that large fluxes of protons (ions) could 
also be extracted from the neutron star and accelerated in the parallel electric fleld in the 
magnetosphere. The initial Poynting flux can be converted into particle thermal and kinetic 
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energy well within the termination radius. Since both electrons and protons are basically 
accelerated by the low frequency electromagnetic wave generated by the pulsar, they will 
be accelerated to the same relativistic speed as they are bound by the strong electrostatic 
force. Hence, the Lorentz factor of the electrons and protons are the same, leading to the 
result that the protons may have carried away most of the spin down power (Coroniti 1990). 
Thus, we can obtain a form for the spin down power from 

E ~ 7V7^mpC^ (11) 

where N is the outflow current from the surface. This current should be of the order of the 
Goldreich-Juhan current (Goldreich & Juhan 1969) given as TV ~ 1.35 x lO^^EigP^^s"^. 
Since E ~ lO^^Sfs^^^erg s-^ we find N oc leading to 7^ oc E^/'^. Therefore, we 

obtain the relation oc E'^/^, where p generally varies between 2 and 3. The relation 
deduced from the non-pulsed X-ray luminosity and spin down power of observed pulsars 
may result from a relatively high electron energy index in the nebula. 



3.2. Statistical properties of pulsed and non-pulsed spectra 

The spectral observations of the X-ray emission provide additional evidence in support 
of the pulsar wind nebula interpretation. In particular, the photon indices of both the 
pulsed and non-pulsed emission components, where determined, for X-ray pulsars are listed 
in Table 2. In this data compilation, the observational results are taken from the ASCA data 
and also from the recent observations by BeppoSAX, Chandra and XMM-Newton. In Figure 
4, the photon index of the pulsed component is plotted versus the ratio of the isotropic 
pulsed X-ray luminosity to pulsar spin down power, Lp^\/E. The results reveal that the 
photon index of the pulsed emission component generally varies from 1.1 - 1.9 (except for 
PSR 0218+4232). Such a range is consistent with the pulsed component originating in 
the pulsar's magnetosphere (Wang et al. 1998; Cheng, Gil, & Zhang 1998). The values of 
Lpui/E vary over a wide range from 10~^ to 10~^ with no apparent correlation of the photon 
index with respect to efficiency. We note that the actual pulsed X-ray luminosities depend 
on the beaming of the pulse and, therefore, the efficiency of spin down power to pulsed 
X-ray luminosity is uncertain. The real efficiency may only scatter around an apparent 
efficiency, which is based on an assumed constant solid angle. Unless the beaming conspires 
with the photon index to produce a correlation, the rough distribution of the data presented 
in Figure 4 would not be expected to change. In summary, the efficiency and the photon 
index do not appear to show any correlation, which is consistent with the model predictions 
of Cheng & Zhang (1999). On the other hand, Cheng & Zhang (1999) suggest that the 
photon index sensitively depends on the local properties of the stellar magnetic field. 
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We also illustrate the relation between the efficiency of pulsar spin down power to 
non-pulscd X-ray luminosity, L^p^i/E, and the photon index of the non-pulsed component 
in Figure 5. It can be seen that the photon index of the non-pulsed emission component 
shows two distinct regimes. In particular, there is a grouping of data points with Fnpui ^ 2 
for low efficiencies 10~^) and a separate grouping of data points with Fnpui 2 at higher 
efficiencies. This tendency may reflect a larger photon index in the fast cooling regime in 
comparison to the slow coohng regime as discussed in §2.3. Specifically, we showed that the 
fast cooling regime and slow cooling regime is separated by {pc/^x) < 1 or > 1 where Uc is 
the cooling frequency and vx is the observed X-ray frequency respectively. In the former 
regime, the photon indices are in the range from 2-2.5 and in the latter regime, the photon 
indices range from 1.5 to 2. Typically, the slow cooling region has a lower efficiency than 
the fast cooling region by a factor {uc/^xY^'^ « 1- The actual efficiency of an individual 
pulsar depends on its corresponding parameters (e.g. p, 7^,i?,ee). Evidence in favor of 
such a relation is suggested by the difference between the Crab, which is in the fast cooHng 
regime with an efficiency of ~ 0.03 and Geminga, which is in the slow cooling regime with 
an efficiency of ~ 10~^. 

3.3. The termination radius and non-pulsed luminosity 

In Table 3, the termination radii derived from the images of some identified pulsar 
wind nebulae are listed. Although the determination of the radius is imprecise, we have 
adopted two procedures for its estimation. For the nebulae which have ring structures (e.g.. 
Crab nebula), the termination radius is chosen as the size of the inner ring. On the other 
hand, for the extended sources with a bright point-like central sources, the termination 
radius is taken to be the size corresponding to the scale in which about 90% of the observed 
counts are included. For comparison, the theoretical termination radii for some nebulae are 
estimated according to Eqs. (7) (for the Crab nebula) and (9) (other nebulae) as discussed 
in §2.2 and also listed in Table 3. 

With the knowledge of the termination radius of the nebulae, the X-ray luminosity 

can be estimated assuming typical values for the parameters (7^, p, eg, e^). The observed 
and estimated X-ray luminosities of these compact nebulae are given in Table 3. For these 
estimates, we have taken /?s as the observed values, and assumed 7^ ~ 10®, p = 2.2, ~ 0.5, 
and ~ 0.01 for all pulsars. We flnd that the model luminosities and the observed 
luminosities are in rough agreement with each other to within a factor of 4, which taking 
into account the possible variations of these parameters amongst the pulsars, suggests that 
the simple one-zone model provides a reasonable approximation for the X-ray luminosity 
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level. 



4. Normal Pulscirs and X-ray Tails 

The X-ray emitting region with a characteristic frequency (vx = ^2m^ ) exhibit a 
tail-hke spatial structure provided that the pulsar velocity exceeds that of the termination 
shock front and the nebula magnetic field is sufficiently low. In this case, the distance 
traversed by the pulsar within the synchrotron cooling timescale can be taken as a lower 
limit of the elongation length. Specifically, the synchrotron cooling time in the X-ray band 
is Tc = 67rmec/7(7T-B^ ~ 10^i?^Q^(/iz/x/keV)~^/^ s where B^a is the magnetic field in the 
emission region in milligauss. Thus, the typical cooling time is ~ 10^^ s for Bmc = 0.01, 
and the length of the X-ray elongated feature is about / ~ VpTc ~ 10^^ cm for a pulsar 
moving at a velocity of 100 km s~^ with respect to the interstellar medium. 

The X-ray images of some pulsar wind nebulae, indeed, reveal extensive X-ray tails. 
In Tabic 3, several sources with X-ray tail features are included. Seven cases arc thought 
to be the consequence of a bow shock formed by the high velocity pulsar (e.g., Gcminga, 
PSR B1823-13, B1757-24, B1957-20, J1747-2958, J1124-5916 and B1853+01). The two 
pulsars J1930-I-1852 and B0453-685 exhibit elongated structures and may also be bow shock 
structures. On the other hand, the Crab and Vela nebulae exhibit an outflow structure 
(X-ray jet). 

4.1. Geminga 

Geminga is one of the most studied pulsars because it is one of the very nearby pulsars 
and is the first pulsar detected in optical. X-ray and gamma-ray, but not in radio (Bignami 
& Caraveo 1996) wavelengths. Based on ASCA observations, Saito (1997a) found that 
the non-pulsed X-ray luminosity was ~ 3.6 x lO^^erg s~^ and the pulsed X-ray luminosity 
was ~ 4.6 X lO^^erg s~^ in the 2-10 keV band. Based on more recent XMM-Newton 
observations, Caraveo ct al. (2004) has reported the detection of a power law component 
(~ 7.7 X 10^*^ erg s^^ in 2-8 keV band), which consists of the pulsed as well as the non-pulsed 
component. The pulsed component in 2-6 keV is ~ (33.2 ± 4.5)%, which is consistent with 
ASCA observation. However, it is still not clear if the non-pulsed component originates 
from the pulsar or the nebula. The recent Chandra observations of the Geminga pulsar 
show a X-ray bow shock structure(Caraveo et al. 2003). The image of the field including 
the Geminga pulsar reveals a compact point-like X-ray source of size 20" or (6 x 10^^ cm) 
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with two elongated tails extending to 2' (Caraveo et al. 2003). Since the X-ray tails have 
a non-thermal X-ray spcctnim (F ~ 1.6), these features may be produced by synchrotron 
emission in the bow shock between the pulsar wind and the interstellar medium. The 
spin down power of the Geminga pulsar is relatively low compared to other pulsars, with 
E ^ 3.2 X lO^^erg s~^. Assuming a proper motion velocity of Vp ~120 km s~^ (Bignami 
& Caraveo 1993) and a number density of 1 particle cm~^ in the interstellar medium, the 
termination shock radius of the pulsar wind nebula is ~ 4 x 10^^ cm. This scale is consistent 
with the observational constraint on the compact X-ray nebula for an assumed distance of 
~ 160 pc. The point-hke source has a X-ray luminosity of 1.2 x lO^^ergs s~^ in the energy 
range from 0.3 to 5 keV, and the total luminosity from the two tails is ~ 6.5 x 10^^ ergs 
s~^ (Caraveo ct al. 2003). The pulsed and non-pulsed emission of the point-like source was 
not resolved in the XMM-Newton observations, but the results from §3 (see Table 1), reveal 
that the non-pulsed luminosity of Geminga is 3.6 x lO^^erg s^^. Interpreting this luminosity 
as arising from the compact X-ray nebula, we find that Lneb,x/-E' 10^^ which is low 
compared to other pulsar wind nebula. We suggest that the nebula is in the slow cooling 
regime where the photon spectral index, T — {p + l)/2. The observed photon index of the 
tail is 1.6, corresponding to p = 2.2, and its observed length (~ 4 x 10^^ cm) suggests that 
the magnetic field of the surrounding medium is weak {B ~ 1.5 x 10~^ G, see also Caraveo 
et al. 2003). Thus, the inferred fractional magnetic energy density is ~ 10~^. Assuming 
a fractional electron energy density of eg ~ 0.5 yields a nebula luminosity of ~ 10^°erg s~^ 
in the 2-10 keV band, which modestly overestimates the contributions from the compact 
nebula and X-ray tails by a factor of ~ 2.5. 

4.2. PSR B1757-24 

Another example where the model can be compared to the observational data is the 
X-ray tail associated with PSR B1757-24 observed by Kaspi et al. (2001) with the Chandra 
satellite. PSR B1757-24 is a 124 ms radio pulsar discovered near the supernova remnant 
(SNR) G5.4-1.2 (Manchester et al. 1985) with a spin down power ~ 2.6 x lO^^erg s~^ 
(Manchester et al. 1991). Assuming an association between the pulsar and G5.4-1.2, and 
that the pulsar's characteristic age is a good estimate of its true age, a proper motion of 
75 mas yr~^ is inferred. This corresponds to a transverse space velocity of Vp ~ 1800 km 
s~^ for a distance of 5 kpc (Frail, Kassim, & Weiler 1994). However, recent interferometric 
observations have failed to detect the imphed proper motion (Gaensler & Frail 2000), which 
suggests that the pulsar is older than its characteristic age or that the assumed pulsar birth 
place is incorrect. Here, we take the proper motion velocity of PSR B1757-24 to be Vp ~ 590 
km s~^ (Gaensler & Frail 2000). The observed 2-10 keV luminosity is about 2 x lO^^erg s~^ 
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corresponding to a conversion efficiency of 10~^, witfi tfie X-ray tail extending to ~ 20" (0.5 
pc for tfie distance of 5 kpc). The observed plioton index, F ~ 1.6, yielding p = 2.2 in the 
slow cooling regime. Assuming an interstellar medium number density, n ~ 1 cm~^, the 
termination shock radius is Rg ~ 10^'' cm, with the known velocity Vp ~ 590 km s^^. Taking 
the typical pulsar wind nebula parameters, eg ~ 0.5, 63 ~ 0.01 (corresponding to ~ 20/iG 
at Rg), and 7^ ~ 10^, we find the X-ray luminosity from 2-10 keV is ~ 3 x 10^"^ erg s"^ using 
the one-zone model. The predicted X-ray tail length is given by / ~ Vptc ~ 2 x 10^^ cm. 
Both these estimates for the X-ray luminosity and tail length are in approximate accord 
with the observed values. 



5. SUMMARY 

The non-thermal non-pulsed X-ray emission of rotation powered pulsars has been 
investigated in the context of emission from a pulsar wind nebula. We have confirmed that 
this emission can significantly contribute to the total X-ray emission, thereby, increasing 
the X-ray emission above that produced in the pulsar magnetosphere. A reexamination 
of ASCA data for pulsars separated into the non-pulsed and pulsed emission components 
indicates that their non-thermal components, individually, are correlated with spin down 
power. In particular, the dependence is steeper for the non-pulsed component with 
Lx oc i^i 'iio i as compared to the relation, oc £'i-2±0-08 f^^. ^j-^^ pulsed component. Our 
results are similar to that discovered by Becker & Triimper(1997), who found that the 
non-thermal (pulsed and non-pulsed) component satisfied the linear relation, cx E, 
based on ROSAT data in 0.1-2.4 keV range. In this paper, we have used the ASCA data 
in the 2-10 keV range, which includes the contribution of the additional thermal pulsed 
component with /cT ~ 1 keV from the polar cap (Cheng & Zhang 1999). In this case, the 
dependence of the X-ray luminosity on spin down power is steeper. Because of the poor 
angular resolution of ASCA, the possible contribution of non-thermal non-pulsed emission 
from the magnetosphere is mixed in with that of the nebula. 

Within the framework of a one zone model for the pulsar wind developed by Chevalier 

(2000), the power law relation between the non-pulsed luminosity and spin down power 
provides an important diagnostic for the electron energy distribution in the shock resulting 
from the interaction of the relativistic wind with the interstellar medium. Specifically, 
Lx oc E^^'^ where p is the power law index of the electron energy distribution, suggesting 
that the index characterizing these pulsars is relatively high, p ~ 2 — 3. The observational 
results also show that the conversion efficiency of spin down power to X-ray luminosity is a 
function of spin down power varying from 10~^ to 0.1. The efficiencies in the lower (upper) 
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end of the range can be understood in terms of emission in the slow (fast) coohng regime 
where the X-ray frequency is less (greater) than the electron synchrotron cooling frequency 

The spectra from pulsar wind nebulae are distinctly non-thermal, described by a 
photon power law index, T — (p + 2)/2 or T — (p + l)/2 for the fast or slow cooling 
regime respectively. This implies a correlation between the photon index and efficiency of 
conversion of spin down power to non-pulsed X-ray luminosity with steeper photon power 
law indices associated with higher conversion efficiencies. Although the photon power law 
indices inferred from the pulsars studied have large uncertainties, there are hints from the 
observed data that such a correlation exists. We note that the low values of F and L^/E 
inferred for the extended emission observed from Geminga and B 1757-24 are consistent 
with such an interpretation. Future observed determinations of the photon power law 
indices could further verify this correlation, thereby, providing an additional diagnostic tool 
for probing the pulsar and its environment. 

Finally, the pulsar wind nebula interpretation may also apply to the spatially extended 
X-ray emission observed in binary systems containing a pulsar. For example, the extended 
X-ray emission associated with the millisecond pulsar B1957-I-20 (Stappers et al. 2003) may 
be interpreted within such a framework. The application of such models to recycled pulsars 
in the Galaxy will be the subject of a future investigation. 

We are grateful to the anonymous referee for the fruitful suggestions. This work is 
partially supported by the NSF through grant AST-0200876, by a RGC grant of the Hong 
Kong Government, and by the National Natural Science Foundation of China under grant 
10273011. 
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Fig. 1. — The total X-ray luminosity (2-10 keV) from ASCA observations versus spin-down 
power of 23 X-ray pulsars. The solid line is Lx = lO^^^-E^^^ as used by Saito (1998), 
and the dashed line represents Lx = 10~^E as used by Becker & Triimper (1997) to 
summarize ROSAT observations of central objects and pulsars. The best fitting function 
is Lx oc £'i-35±o.2^ where the error in the power law exponent represents ilcr corresponding 
to the scatter in the observed data points. 
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Fig. 2. — The pulsed X-ray luminosity (2-10 keV) from ASCA observations versus spin-down 
power of 23 X-ray pulsars. The solid hne is Lx = 10~^^£^^/^, and the dashed hne represents 
Lx — 10~^£'. The relation between the pulsed component and spin-down power cannot be 
described by both the two formulae. The best fitting function is shown as the dotted line, 
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Fig. 3. — The non-pulsed X-ray luminosity (2-10 keV) from ASCA observations versus spin- 
down power of 23 X-ray pulsars. The solid hne is Lx — 10~^^£'^/^, and the dashed hne 
represents Lx — 10~^£'. The best fitting function is Lx oc 
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Fig. 4. — The photon index of the pulsed emission component of spin powered X-ray pulsars 
versus 77pui = -^pui/ the ratio of the isotropic pulsed X-ray luminosity to spin down power. 
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Fig. 5. — The photon index of the non-pulsed emission component of spin powered X- 
ray pulsars versus 77npui = L^pui/E, the efficiency of spin down power to non-pulsed X-ray 
luminosity. The dashed line corresponding to F = 2, which approximately separates the fast 
cooling regime from the slow cooling regime, is shown for convenience. 
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Table 1: Characteristics of spin powered pulsars and their X-ray luminosities observed by 
ASCA 
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Note. — The first column PSR is the pulsar name, P is the spin period, P is the period derivative, d is 
the distance of the pulsar from us in units of kpc. The luminosity is in units of erg s~^. E is the pulsar's 
spin-down power. Lx,tot is the total X-ray luminosity observed by ASCA (2-10 keV), -C/x.pui is just the pulsed 
X-ray luminosity, -^^x.npui is the non-pulsed luminosity. Because the X-ray luminosities are all calibrated by 
ASCA, we can assume a 10% uncertainty for the total flux measurement, and except for B1706-44, B1046-58, 
J1610-50 and J0437-47 which have a 50% uncertainty for the pulsed fraction, the others have 30% fraction 
uncertainties. 

References: 1. Torii et al. 2001; 2. Torii et al. 1997; 3. Torii et al. 1999; 4. Saito et al. 1997a; 5. Helfand, 
Gotthelf & Halpern 2001; 6. Pavlov et al. 2001; 7. Halpern & Wang 1997; 8. Cavaveo et al. 2003; 9. 
Pivovaroff ct al. 2000; 10. Wang & Halpern 1997; 11. Hirayama et al. 2002; 12. Shibata ct al. 1997; 13. 
Halpern et al. 2001; 14. Marshall et al. 1998; 15. Wang & Gotthelf 1998; 16. Gotthelf et al. 2000; 17. 
Takahashi et al. 2001; 18. Sakurai et al. 2001; 19. Saito et al. 1997b. 
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Table 2: X-ray pulsed and non-pulsed component properties of spin-powered pulsars 
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Note. — The spin down power is in units of erg s"^. 77pui = Lpu\/E and rynpui = inpui/S arc the X-ray 
efRciencies of spin down power to pulsed and non-pulsed X-ray luminosities separately. Fpui and Fnpui are 
the pulsed and non-pulsed photon indices respectively. 

References: 1. Gacnslcr et al. 2003; 2. MarshaU et al. 1998; 3. Kaaret et al. 2001; 4. Marsdcn ct al. 1997, 
5. Mineo et al. 2000; 6. Torii et al. 2001; 7. Torii et al. 1999; 8. Willingale et al. 2001; 9. Pavlov et al. 
2001; 10. Halpern & Wang 1997; 11. Cavaveo et al. 2003; 12. Shibata et al. 1997; 13. Takahashi et al. 
2001; 14. Saito et al. 1997b; 15. Halpern ct al. 2001. 16. Gotthelf & Kaspi 1998; 17. Finley et al. 1998; 
18. Kaspi et al. 2001; 19. Murray ct al. 2001; 20. Wang & Halpern 1997; 21. Greiveldinger et al. 1996; 22. 
Hessels et al. 2004; 23. Gaensler et al. 2004; 24. Petre et al. 2002. 
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Table 3: X-ray properties of pulsar wind nebulae 
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Note. — The luminosity is in units of erg s~^. and Rf^ are the observed and predicted termination 
radius of the pulsar wind nebulae. i°pui and are the observed and theoretical X-ray luminosities of 

pulsar wind nebulae, n is the number density of the medium around the pulsar. Vp is the pulsar proper 
motion velocity, I is the length of the observed X-ray elongated structure. 

^ Seven fast moving pulsars in the bow shock show the X-ray tail features. I of the pulsars Crab and Vela 
corresponds to the outflow scales of the nebulae. For B0453-685 and J1930-I-1852, I how the scales of the 
elongated structure. 
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1994; 5. Weisskopf et al. 2000; 6. Pavlov et al. 2001; 7. Cavaveo et al. 2003; 8. Halpern et al. 2001. 9. 

Gotthelf & Kaspi 1998; 10. Finlcy et al. 1998; 11. Kaspi et al. 2001; 12. Murray ct al. 2001; 13. Stappers 
et al. 2003; 14. Hessels et al. 2004; 15. Gaensler et al. 2004; 16. Camilo et al. 2002; 17. Hughes et al. 2001; 
18. Petre et al. 2002; 19. Lu et al 2002; 20. Gaensler et al. 2003b; 21. Romeni & Ng 2003. 



